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Abstract 

> 

, We establish a new approach to calculating spectral statistics in disordered 

conductors, by considering how energy levels move in response to changes in 

the impurity potential. We use this fictitious dynamics to calculate the spee- 
ch i 

tral form factor in two ways. First, describing the dynamics using a Fokker- 
Planck equation, we make a physically motivated decoupling, obtaining the 
q | spectral correlations in terms of the quantum return probability. Second, 



X 



from an identity which we derive between two- and three-particle correla- 



tion functions, we make a mathematically controlled decoupling to obtain the 

same result. We also calculate weak localization corrections to this result, and 

show for two dimensional systems (which are of most interest) that corrections 

vanish to three-loop order. 
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I. INTRODUCTION 



Numerous properties of quantum systems can be described in terms of their energy 
spectra. For complex systems an exact determination of energy levels is not feasible, and 
a statistical description becomes necessary. It turns out that the Wigner-Dyson (WD) 
statistics of eigenvalues of random Hermitian matrices describes energy levels in a wide 
variety of different systems [[|. The joint distribution of eigenvalues is dominated by level 
repulsion and is universal in the sense that level correlations depend only upon the symmetry 
of the Hamiltonian while all specific properties of the system are absorbed into the mean 
level spacing, A. A very important feature of WD statistics is that - by construction of 
invariant ensembles of random matrices - spectral properties are independent of eigenstate 
correlations. In real systems such an independence can at best be approximate. It holds, 
however, in the ergodic regime where the entire phase space of a system is explored. If 
a non-ergodic regime is of interest, not only is WD statistics inapplicable but the whole 
concept of the independence of spectral and eigenstate correlations should be re-examined. 

Disordered mesoscopic conductors present a natural ensemble for a statistical description 
- the ensemble of impurity configurations. In this case spectral statistics in the non-ergodic 
regime are very important for both transport and thermodynamic properties of electrons. 
Different regimes in disordered conductors are determined by the energy or time scale, as 
shown in Fig. 1. The ergodic regime involves energy level separations e ^ E c = h/t erg where 
E c is called the Thouless energy, and t erg ~ D/L 2 is the time required for the electronic 
diffusive motion, with diffusion constant D, to fill all phase space, in a sample of size L. 
The quantum limit of this regime corresponds to smaller energy separations, e ^$ A, where 
A is the mean level spacing, and to longer times t >t H , where t H = h/A is the Heisenberg 
time. Note that the ratio E c /A is proportional to the dimensionless conductance g (i.e. 
the conductance measured in the units of e 2 /h), and is large in the metallic phase. The 
diffusive regime involves energies h/t eTg < e ^ fo/t e i, where t e \ is mean elastic-scattering 
time. The largest energies, h/t e \ < e ^ e F , (where e F is the Fermi energy) and the shortest 
times, t < £ e i, correspond to the ballistic regime in which multiple scattering by impurities 
is irrelevant. 

It was first conjectured by Gor'kov and Eliashberg [f|] and then shown by Efetov || 
that spectral correlations in the ergodic regime are described by the random matrix theory 
(RMT) of Wigner and Dyson. Correlations in the non-ergodic diffusive regime which are 
important, in particular, for universal conductance fluctuations were analysed by Altshuler 
and Shklovskii || at leading order in diagrammatic perturbation theory. Their results were 
later reproduced by Argaman et al 0] within the semiclassical approach, using the diagonal 
approximation. 
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FIG. 1. Regimes of energy and time in a disordered metal; here s e i = h/tciA and s F = s F /A. 

We have recently developed || an alternative approach to level statistics in the non- 
ergodic regime which takes into account the inevitable coupling between eigenvalue and 
eigenstate correlations, and can be extended beyond the region of validity of the pertur- 
bative technique. Our approach is based on the idea of parametric motion through the 
ensemble of disordered Hamiltonians: we treat the energy eigenvalues as particles with fic- 
titious dynamics induced by changing some parameter of the Hamiltonian. This dynamics 
takes the form of Brownian motion in a fictitious time r related to the parameter being 
changed. Originally, this idea was employed by Dyson H] in the context of RMT. Later, 
in the context of the semiclassical description, Pechukas |10[ used motion along a smooth 
path in the space of Hamiltonians to generate fictitious dynamics of a different kind. Both 
Dyson and Pechukas were interested in level dynamics primarily as a way of generating the 
level distribution. By contrast, a number of recent authors notably Szafer, Simons 



and Altshuler [|X2|JX3[1 , have investigated the dynamical problem in its own right, calculating 
parametric statistics: eigenvalue correlations as a function of position in the space of Hamil- 
tonians. In distinction to our approach, eigenfunction correlations have played no role in 



previous work (TT|-|T3|,|T5|,|16[ , an assumption justified only for the ergodic regime. 

We have found || that a treatment based on Brownian motion through the ensemble of 
Hamiltonians provides a unified description of all regimes in disordered conductors, except 
the quantum limit (t > t H , or E < A), which is also beyond the scope of diagrammatic and 
semiclassical approaches. The main result is a new relation explicitly linking the spectral 
correlation function to the quantum return probability for an expanding wavepacket which, 
in turn, is related to a certain eigenfuction correlator. The derivation given in Ref M has 
a limitation: when obtaining a closed Langevin equation describing the Brownian motion, 
we make use of an uncontrolled, although physically transparent assumption. Thus, within 
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the framework of that calculation it is not possible to establish exactly the region of validity 
and accuracy. 

In this paper, we re-derive the relation between spectral and eigenstate correlation func- 
tions using a more explicit procedure: we decouple a certain exact relation between two- and 
three-level correlation functions using the Kirkwood superposition approximation. Then we 
examine the accuracy of this decoupling using perturbative diagrammatic techniques. Re- 
markably, it not only reproduces the results of the diagonal approximation ||[7]], but holds 
well beyond it. We show this to third order in a perturbative expansion in for two- 
dimensional systems with or without time-reversal symmetry. We are therefore encouraged 
to believe that the connection between spectral and eigenstate correlation functions should 
be useful rather generally, and especially for problems where the usual diagrammatic tech- 
nique cannot straightforwardly be applied, such as spectral statistics in the critical regime 
near the Anderson transition. 



II. DEFINITIONS AND MAIN RESULTS 

A convenient way to consider spectral correlations is to introduce fictitious level dynamics 
in response to changing some parameter A of the Hamiltonian. Thus we parameterize the 
ensemble of Hamiltonians as follows 

H(X) = H + XW(r) . (2.1) 

Here both H and H(X) belong to the same symmetry class, and the point A = corresponds 
to some arbitrary choice of one of the many members of the same ensemble. We will specify 
the choice of H and W in the next section. 

We consider in this paper the two-level correlation function (TLCF) and its Fourier 
transform, the spectral form factor, in a disordered conductor described by the Hamiltonian 
Eq. (|2.1| ). Let E n (X) be the energy levels of H(X). We introduce the density of states per 
unit volume (DoS) as 

p(E,X) = j- d J25(E-E n (X)). (2.2) 

The mean level spacing, A, is then related to the mean DoS, p = (p), by A = (pL d ) _1 . The 
TLCF is defined as 

R(s, A) = p- 2 (p(E + sA, X)p(E, 0)) - 1 , (2.3) 

where u = sA is the energy difference between two levels. The mean DoS is practically a 
constant in the entire energy region of interest (as it changes only at scale of order e F while 
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we consider energy windows centered at e F of width not exceeding %/t e \ <C e F ). We consider 
only values of A small enough so that the statistical regime does not change and neither does 
the mean DoS (for large enough L this nevertheless allows arbitrarily large A on the scale 
relevant for parametric correlations). Because of this the TLCF cannot depend on either 
E, or on the choice of the point Hq in the ensemble ( |2.1|) . We define the (dimensionless) 
spectral form factor as 

/oo 
e~ lst/t H R(s, X) ds . (2.4) 
-oo 

Our main result relates the spectral form factor to the quantum return probability p(t) 
of a diffusing electron as follows: 



= (2.M-WrW (25) 
i + (*h P )-i j$ P (f)dt> 



We have obtained this expression for times shorter than the Heisenberg time t H = h/A. 
Here we define p(t) as the probability density for the wave packet, originally created in a 
small volume Vq ~ £ d , to remain in this volume at the time t (£ is the elastic mean free path 
which is a natural coarse-graining size for the disordered metal; however, we could choose Vq 
arbitrarily, provided that £ d > Vq ^ A^ where A F is the Fermi wavelength). The ensemble- 
averaged return probability is related, as we will show later, to the following wave-function 
correlations: 

p(t) = Jd d r (T |^(r)|> n+i (r)| 2 e- i ^-^)*/^ . (2.6) 

It is important to note that, by definition of the wave packet above, the summation here is 
limited to the number of levels TV ~ L d /Vq with energies lying within the energy band of 
width E ~ 1/pVo- 

Equation (|2.5| ) relates the spectral and wave-function correlations. Let us analyse this 
relation in the metallic phase. In the diffusive regime, t e i ~ ^ ~ 4r g , the quantum return 
probability p(t) reduces at leading order to the classical return probability for random walks, 
multiplied by a symmetry factor 2/f3, where (3 — 1, 2 or 4 is the usual index corresponding 
to the orthogonal, unitary, and symplectic symmetry ensembles, respectively [[|: 

^^wktfF-- (2 ' 7) 

Noting that in the ballistic regime, t < t e \, p(t) saturates at po(t e {) ~ l/£ d , one sees that 
the integral in the denominator of Eq. fl2.5p is of order (t e \A/h)(L/£) d ~ g^ 1 for d > 2, 
and of order g^ 1 hi(t/t e \) for d = 2. It is well known that such an integral describes a 



weak localization correction to conductance and other physical quantities JT^|T^]. On the 
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other hand, the quantum return probability p(t) contains weak localization corrections itself. 
Neglecting these corrections in both the numerator and denominator of Eq. ( |2.5| ), we reduce 
it to 



K (t) = (2irhp)- 1 \t\p (t). (2.8) 

To leading order, this expression is also valid in the ergodic regime, t C r g ~ t "C t H , where the 
classical return probability saturates at (2/ (3)L d so that the second term in the denominator 
in Eq. ( |2.5|) is of order t/t H <C 1 and may be neglected. We should not expect Eq. ( |2.5|) to 
be correct in the quantum limit, t 3> t H , as we have derived it under the assumption that 
the opposite inequality holds, as will be seen later. Indeed, in this regime Eq. ( |2.5| ) gives 
the saturation of K(t) at 1/2 instead of the correct limiting value K(t) = 1. 

Equation ( |2.8[ ) coincides with the result obtained by Argaman et al 0, using the di- 
agonal approximation in semiclassical periodic-orbit theory. The Fourier transform of this 
expression corresponds to the TLCF obtained originally by Altshuler and Shklovskii ||. In 
the diffusive regime, R(s,0) ~ A d g- d / 2 s d / 2 - 2 , where A d is a numerical coefficient which is 



zero for d = 2 [19], and in the ergodic regime, R(s, 0) ~ —1/s [pO 



The second relation we obtain between K(t) and p(t) is: 

K(t) + (rrhp)- 1 / K(t - t')p(t')dt' = (2nhp)- 1 \t\p(t) (2.9) 



CH 

which we can see is very similar to Eq. ( |2.5| ). The latter is obtained from a diagrammatic 
analysis discussed in section (VI). The new feature is that we have a convolution of K(t) 
and p(t) which occurs because the decoupling is in cu-space rather than t-space. In 2d, up 
to 3-loop order in perturbation theory, both these relations correctly reproduce the TLCF. 
It seems to us quite remarkable that a relation derived from a phenomenological model of 
energy level dynamics could be exact to 3-loop order. We note that the 2d case is expected 
to be a good model for the behaviour of the system for d > 2 at the mobility edge. From the 
point of view of a power-counting analysis of the properties of K(t) at the mobility edge, 
both Eq. fl2.5|) and Eq. fl2.9|) should work equally well. 



III. RANDOM WALKS THROUGH THE ENSEMBLE OF HAMILTONIANS 

We have established elsewhere || the relation ( |2.5| ), using a Langevin equation to describe 
the motion of levels on the energy axis in response to a random walk through the the ensemble 
of Hamiltonians. In contrast to eigenvalue statistics in RMT where the Brownian motion 
ideas were originally applied [^,|T^], the Langevin equation for level motion is not closed, and 
certain assumptions are required to solve it. In the following section, we will re-derive Eq. 
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Q2.5|) with the help of a different approach based on the decoupling of a certain exact relation 
between two- and three- level correlation functions. Before doing this, however, it is useful 
to analyze the Brownian motion picture within the Fokker-Planck scheme. Although the 
Langevin and Fokker-Planck schemes are in principle equivalent, the assumptions required 
in order to make the description closed are different. The Langevin scheme of Ref. |J is 
physically more transparent. The advantage of the Fokker-Planck scheme which we develop 
here is that the approximation made are more closely related to those analysed in subsequent 
sections. 

We consider paths of two types through the ensemble of Hamiltonians Q2.1|) which, for 
free electrons in a random potential, have the form 

h 2 

H(X) = -— V 2 + U(r) + XW{r) . (3.1) 
2m 

Here both U(r) and W(r) are chosen to be of Gaussian white-noise form with zero average 
and 

(U(r)U(r')) = -±—6(r-r>), 

2npt el (3.2) 

(W(r)W(r')) = v 2 L d 5(r-r') . 

The first type of path is a straight line through the ensemble, and generated by varying A in 
Eq. The second type is a Brownian path through the ensemble, parameterized by the 
fictitious time r, generated in the following way. 

H(t) = H + f T dr' V(t', r) . (3.3) 
Jo 

We take V(r, r) to be Gaussian distributed with zero average and 

(V(t, r)V(r', r'))=v 2 L d 5(r - r') 5(r-r') . (3.4) 

Referring to Fig 2, one sees that the two ways of exploring the ensemble are equivalent if 
one makes the identification r = A 2 . 

In our derivation we will use averages over both W or, equivalently, V and then over H Q , 
and we must here discuss the role of each. The average over all possible perturbations, W, 
will be necessary to derive the equation of motion for the of the energy levels, and thence the 
density of states. We can then obtain correlation functions for energy levels of the system at 
different parameter values by averaging over the starting point Hq. Such functions should 
then depend only upon energy and parameter differences by homogeneity. 

The first step in our derivation is to obtain the equation of motion for the joint probability 
density function (JPDF) of energy levels, P({E n }, r). We use perturbation theory to second 
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FIG. 2. Smooth [H (A)] and Brownian [H(t)] paths through the space of Hamiltonians 

order to calculate the change of E n (r) in response to the evolution from r to r + 5t. After 
averaging over W we obtain 

(f)V)l ' 2 £ p M (3 ' 5a) 

m _, n E n (r) - E m {T) 

[SE n (r)SE m (r)) = v 2 c nm (r) (3.5b) 
where 

c nm {r)=L d fd d r\Mr,r)\ 2 \^ m (r,r)\ 2 , (3.6) 



5EJr) 



and vp n (T,r) are the corresponding eigenfunctions of H(t). Before we can use the above 
equations to derive a Fokker-Planck equation for the JPDF, P({E n },r), we must make an 
assumption. We replace c nm {r) by its average over the ensemble of H , This amounts to 
ignoring correlations between eigenvalues and eigenvectors. We take the disorder average to 
be a function only of the energy difference, u = E n — E m , within the window of interest: 

(cnm(r)) = c(u) . (3.7) 

Furthermore, we assume that in the Fourier transform of the wavefuntion correlator c(u) we 
may neglect the correlations between the eigenvectors and eigenvalues so that 



duj 
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to the return probability of a diffusing electron, Eq. (|2.6|). To this end, consider a wavepacket 
made from the eigenstates of H{r) and concentrated initially in a volume Vq of order l d near 
the origin (since r plays no role, we suppress it as a label in the following): 



^(r,t) = Aj2M0TMr)e 



-iE n t/h 



Here the summation is limited to M ~ {L/£) d levels with energies \E n \ < l/pl d , and the 
normalization constant is A 2 = L d /N '. The ensemble-averaged return probability p(t) = 
(|^(0, t) | 2 ) is given by 



pit) = A 2 ^(\4> n (0)\ 2 \^ m (0)\ 2 e- i{E "- Em)t/h ) = {^J2\^n(0)\ 2 \ij n+ i(0)\ 2 e~ i{E "~ En+l)t/n 

nm \ I I 

where we have used the fact that the first sum above depends only on the difference \n — m\. 
Noticing also that the ensemble-averaged quantity is spatially homogeneous, we reduce this 
expression to that given in Eq. ( |2.6| ). Comparing this to the definition of c nm , Eq. ( |3.6| ), we 
obtain for t > 

P(t) = T d (^ Cn,n + ie- i{En - E ™ +l)t,h ^ • (3.9) 
On the face of it, this coincides, up to a constant factor, with C(t), the Fourier transform 



of c n>n+ i, introduced in Eq. (|3.8|) . There is, however, an essential difference: C(t) is defined 
by the Fourier sum containing all the levels (say, up to e F ), while the Fourier sum for p(t) 
contains only the levels within the bandwidth E <C e p . When \E n — E n+i \ > Eq the two 
levels are practically uncorrelated, and c(u) = 1 for u > E , so that C(t) contains a 5-like 
function for t near zero. As we are not interested in an exact description at the ballistic 
time scale, we can represent the relation between C(t) and p(t) as follows: 

pit) = 2nhpC{t) , t > 0; 
r t i i f t (3.10) 

We also note that the definition of H(t) in Eq. ( |3.3| ) causes the energy levels to move 
away from each other indefinitely as parametric time increases. To overcome this problem 
we introduce a rescaling term, —5rU{E ri ) to the r.h.s. of Eq. ( p.5a|) . This U{E n ) can 
be thought of as a Lagrange multiplier, and it will be set later on by the condition that 
correlation functions can depend only on differences in parametric time. 

With these considerations, starting with Eqs. ( |3.5|) we end up with the Fokker-Planck 
equation for JPDF: 

1 OP d ( dU \ ^ d 2 

V 2 ^ ~~\dE~ VdE n ) + h mM: [CnmP) 11 1 1 ] 
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where the drift potential term, U is the sum of one-particle and two-particle potentials, 



U{{E n }) = J2 U(E n ) + f(En - E m ). 



(3.12) 



The one-particle potential arises from the energy rescaling described above, and may be con- 
sidered as a confinement potential for a one-dimensional gas of fictitious particles interacting 
via the two-particle potential, f(u)), which is related to c{u) by 

df(u) c(u) 



du 



(3.13) 



We see that both the drift potential and the diffusion term in the Fokker-Planck equation 
are expressed in terms of the function c(u), which we have shown to be related to the 
return probability p(t). The off-diagonal diffusion terms in Eq. (|3.11|) , which are due to 
eigenfunction correlations, mean that the static solution does not have a simple Gibbs form. 
In fact we cannot write down its solution in closed form at all. The absence of a simple 
static solution to the Fokker-Planck equation fl3.11|) is an important difference between the 
current problem and the Brownian motion approach to RMT || where such a solution yields 
the exact JPDF. However, the JPDF contains much more information than we require; for 
our purposes it is sufficient to study the equation of motion for the density of states, which 
can be written in the form 



where 



p(E,r)=L- d Y,S(E-E n ), (3.14) 

n 

means averaging over the JPDF P({E n },r). Following the procedure of Dyson 



!| and Pastur [22] we obtain, 



1 dp(E, r) d 2 . . . , „ . d 
= ^[c(0)p(£,r) + — 



dr 



dE 



— + L JdE P2 (E,E,r) — 



where 



p 2 (E, E', r) = LT Yl S ( E ~ E n)^ E - E n 



(3.15) 



(3.16) 



Since we are only interested in level correlations in energy windows small compared to the 
scale at which p varies, the first term in Eq. (|3.15|) is negligible. We rewrite the last term in 
Eq. (glgD in terms of the static TLCF, R(E - E'), 



p 2 (E, E', r) = p(E, r)p(E>, r)[l + R(E - E% 



(3.17) 



After substitution into Eq. ( |3.15 ) we see that the integral over ppR is dominated by a region 
of relatively small energy differences since the product Rf falls off rapidly. In this region p is 
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roughly constant, and the integral vanishes by oddness of the integrand. We have therefore 
arrived at the equation 



1 dp(E,r) _ d 
" dE 



dr 



p(E, r) A (u(E) +L d J dE'p(E>, r)f(\E - E'\ 



(3.18) 



This is a non-linear equation, and to proceed further we must linearize it. The static 
solution of Eq. ( |3.18|) is just the equilibrium density of states, p eq (E), and we perform 
our linearization around by expanding around p eq (E), to give the following equation for 
p(E,r) = p(E,r)-p eq (E): 



1 dp(E,r) 1 d 



J-A / dE'f(\E- E'\)J-p(E' ,r) 
A dE J Ml u dE' HK ' ; 



dr 



where we have approximated p eq (E) by p. Multiplying both sides of this equation by p(E", 0) 
and averaging over the starting point H gives us the evolution equation for the TLCF, 
R(E,t), 

d r dJ 



d_ 



(3.19) 



where we have fixed the units of r by setting v 2 = A/nh. Eq. (|3.19|) can then be solved by 
taking the Fourier transform to yield the result for the spectral form factor 



K(t, t) = K(t, 0) exp 



M(t) 



tT 



(3.20) 



where M(t) = 2tf(t). From the definition of f(u), Eq. (|3~T3|) , we see that M(t) is related 
to C(t), the Fourier transform of c(u), by 



M(t) = 2 f C(t')dt' =1+ ' 
Jo 



p(t')dt . 



(3.21) 



nhp Jo+ 

Although Eq. (|3.20|) gives the parametric dependence of K(t, r) in terms of a function 
M(t) related to eigenfunction correlations, we still do not know K(t, 0). To relate K(t, 0) to 
eigenfunction correlations we introduce a Ward identity as follows. Similarly to the TLCF, 
Eq. ( |2.3|) , we define the current-current correlation function: 



C(s, X) = A 2 J2 (EnW E m (0) 5(E + sA - E n (X)) 5{E - E m (0)) 



(3.22) 



The assumption that both the correlation functions depend only upon energy and parameter 
differences leads to the Ward identity 



d 2 C(u, A) d 2 R(uj, A) 



du 2 



d\ 2 



(3.23) 
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and thence to the relation 

dK(t, t) A 



-- JduC(u,0)e~ iwt/n . (3.24) 



dr 

Finally we can relate C(u, 0) to eigenstate correlations assuming, as above, that we can 
ignore higher order correlations between eigenvalues and eigenstates: 

C{u, 0) = A 2 ]T (Wnn W mm 5(u - E n ) S(E m )) « v 2 c(u) (3.25) 



from which it follows that 



K ^^Q 2 wr (3 - 26) 



With allowance for Eqs. ( 3 . 1 0|) and 3.21 ), this is equivalent to the relation (|2.5|) obtained 



within the Langevin picture of Ref . || . 

The crucial assumptions used in the Brownian motion approach were in the lineariza- 
tion of appropriate equations, and in the neglection of higher order correlations between 
eigenvalues and eigenstates. We believe that these assumptions are reasonable provided 
that one considers only behaviour at energy scales much larger than the mean level spacing. 
However, it is not possible to establish exactly their region of validity and accuracy within 
the Brownian- motion approach developed here and in Ref. |J. This we will do in the next 
section using an alternative approach. 



IV. EXACT RELATIONS BETWEEN THE RETURN PROBABILITY AND 
HIGHER-ORDER CORRELATION FUNCTIONS 

Our aim now is to re-derive Eq. ( [2.5D using some exact relations which involve the return 
probability by making only explicit assumptions whose region of validity can later be verified. 
We again consider the ensemble of Hamiltonians, each describing a particular realization of 
the impurity potential. Now it will be more convenient to use the representation of Eq. ( |2.1| ) 
where the path through the ensemble is a straight line. We define the Fourier transform of 
DoS as follows: 

/oo 
p{E, \)e~ iEt/h dE = Y, e- iE ^ x)t/h . (4.1) 

As all members of the ensemble are statistically equivalent, averaging over realizations should 
give the same result whatever point along the path has been chosen. Thus one should have 

(K(t, \)K(t', A)) = (K(t, 0)K(t', 0)) (4.2) 
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Now we write for small A 

1 



n(t, A) = n(t, o) + \n(t, o) + ^x 2 n(t, o) + o(A 3 ) 



and substitute this expansion into Eq. (|4.2| ) to obtain 

2 (n(t, 0)7t(t', 0)) + (ft(t, 0)K{t', 0)) + (ft(t', 0)7e(t, 0)) = . (4.3) 



We will use this identity to derive exact relations between spectral correlation functions. It 
follows from the definition fl4.ll) that 



n 



(4.4) 



e -iE n t/h 



where 

|2 



(n| W|m ) 

-EL — E r , 



E n = (n\W\n), E n = Y, 

m 

First consider (TZ(t)lZ(t')) . Averaging over W only, we have (E n E n+ i) = v 2 Lr d c n , n +i- 
Hence, averaging also over Ho, we obtain 

(lim(t>)) = ^Sit+t') (^£c n , n+ ,e-^-W/^ = 2 -^ P (t) , (4.5) 

where we have used Eqs. ( |3.9|) and (|3.10| ) to relate the average above to the return probability 
p(t). The crucial assumption in the derivation of Eq. ( f4.5|) was that of homogeneity in energy 
space which means that (c n , n+ i) does not depend on n. We expect this assumption to be 
valid in the whole energy range of interest since the mean density of states is a disorder- 
independent constant in the energy window centered near e F of width Eq < e F . 

Next, consider (H(t)lZ(t')} , initially averaging only over W. Noting that (it(t)lZ(t') ( 

n(t')(1Z(t)) , as K{t') = K(t, 0) does not depend on W } and 
we obtain from Eq. fl4.4j ) after rearranging the terms in the summation 

e -iE„t/h 



« M \ itv 2 j 1 - e <E n -E n+l )t/h it 

>w hL d V [t^o E n -E n+ i h 
tv 2 ^ , /R ft 



Y^ e - lEn t/n r dt « e i(E n -E n+l )t»/K Cnn+i _ 

tf Jo 
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Hence, averaging on H$ as well and using the same assumption as in the derivation of Eq. 
Q4.5D, we have 

(it(t,o)n(t',o)) = -^^5(t+t r ) fdfl^ y^ Cnn+ie ^-En +l )tyn e -i(E n -E n+mW n\ (46) 

\ / hA Jo \L d lm ' I 

On exchanging t and t', one gets the same expression. 

Now, substituting Eqs. ( |4.5|) and Q4.6|) into Eq. ( |4.3|) , we obtain the following exact 
relation, 

tp(t) = Q(t) = J' dt' l\Y. °n,n+l e^ En - En + l)t ' /% e -i{E n -E n+m )t/h\ ^ 

\ lm I 

where p(t) is the return probability, defined by Eq. (|2.6|) . This exact relation is a starting 
point to investigate the nature of the approximations needed to obtain our Brownian motion 
formula, Eq. (|2.5| ). The exact relation involves the higher order correlations, and we need to 
construct some decoupling procedure. This will be done in the next section. Then we will 
use the diagrammatic technique to analyze the accuracy of the decoupling. 



V. DECOUPLING OF THE HIGHER-ORDER CORRELATIONS IN w-SPACE 

AND t-SPACE 



It is seen from Eq. ( p.6|) that pit) involves correlations of two energy levels, whereas, from 
Eq. (|4.7p , Q(t) involves correlations of three energy levels. We should therefore assume that 
Q(t) can be factorized in terms of two- level correlations. In order to see how to factorize 
Eq. ( [4.7|) , let us note that, using homogeneity of the energy space and the definitions of the 
DoS and TLCF, Eqs. Q and (|p|), one can represent the TLCF (at A = 0) as 



R(lu) = R(sA) = A - E n + E n+m )J - 1 . 

Then, from the definition ( ^.4| ) one obtains the following representation for the form factor: 



K{t) = ( E e -^- E "^ h ) - ^5(t) . (5.1) 



The 5 function above cancels that arising from summation over high levels in Eq. fl5.1|) . 
Naturally, K{t) defined by Eq. (|2.4|) as the form factor for the irreducible TLCF, Eq. ([2. 3D, 
is regular at t — 0. It is seen now that the natural factorization of Eq. ([4.7| ) is 



Q(t) « 2 jT dtl^Y. Cn,n + ie^- E ^'^ (j2 ' 



)t/h 



(5.2) 
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E n+ i E n+m 



weights: |%| 2 |¥ n+1 | 2 



FIG. 3. Structure of the three level correlation function Q(t). We see that if in factorization 
we consider the correlation between E n and E n+m , we must also consider that between E n+ i and 



The only feature that deserves comment is the factor of 2 on the r.h.s. of Eq. fl5.2| ). To see 
how this arises we rewrite Eq. ( |4.7| ) using the definition of c n>n+ i, 

Q{t) = fdH Jd d r (jr |^„(r)| 2 \ip n+ i{r)\ 2 e <E n -E n+l )t> /% e i(E n+m - En )t/n\ 

\l,m I 

This can be represented on the energy axis in the schematic form shown in Fig. 3. We 
should allow both for correlations of E n+m with E n and E n+m with E n+[ , since these are 
equivalent. This is made manifest by changing variable in the integral from t' to t" = t — t', 

Q(t) = j*df Jd d r (j2 \Mr)\ 2 |^ n+z ( r )| 2 e ^«+^ s «)*"M e ^ +m -^ + ;)*/^ 

We rewrite Eq. ( |5.2| ), taking into account that (i) Q(t) = tp(t), (ii) the first average in Eq. 
equals 2irhpC(t'), Eq. (|3.8|), and (iii) the second average is equal to K(t), Eq. 

1 tp(t) 

iTihp f*dt'C(t') 

which is, with allowance for Eq. ( p.lOj ), exactly equivalent to the Brownian motion result 
of Eq. (|2.5| ). Now it is clear that the assumptions we have made to derive Eq. (|2.5|) are 
equivalent to neglecting three-level correlations and keeping only two-level correlations. In 
the above derivation, we have disregarded the 5 function coming from Eq. (|5.1|) , as the exact 
relation (|4.7|) has been actually derived from Eqs. ( |4.5|) and (|4.6|) as t 2 p(t) = tQ(t), and this 
5 function enters in combination tS(t). 

For further analysis of the accuracy of Eq. ( |2.5|) , it will be useful to see how the factor- 
ization ( f).2\) arises in the energy representation. We begin by representing Q(t) in Eq. ( |4.7|) 
as follows: 

Q(t) = fdt' fdu f duo' e - iw ' t,/n e lujt/h Q{u', to) , (5.4a) 



m = ^ 1 T^77K (5-3) 



o 



Q(u',u) = L~ d (j2 6 ( E n+i-En-u')5(E n+m -E n -Lj)c n , n+l ) . (5.4b) 
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The function Q(u', u) can now be related to the three-level correlation function Q(E", E', E) 
defined by 



Q(E",E',E)=L- d /Yl HE"-E n )6(E'-E n+l )6(E-E n+m )c nin+l j (5.5a) 

\n,l,m I 

= L d Jd d r (p(E", r)p{E', r)p(E)) = L d Jd d r Q(E", E', E; r) (5.5b) 

where 

p(E',r)=J2\Mr)\ 2 S(E-E n ) (5.6) 

n 

is the local density of states (LDoS). Thus the correlation function involves, by definition, 
correlations of eigenvalues and eigenstates. By the assumption of homogeneity of energy 
space we see that Q(E",E',E) is a function only of energy differences u = E — E" and 
uj' = E' — E", so that we can put E" = without loss of generality. Integrating Eq. (|5.5a|) 
over E" prior and after putting E" to then yields the following identity: 

MQ(uj\uj) = E Q{0,u',u) 

where E is energy bandwidth and M is total number of energy levels. From this identity 
and Eq. ( |5.4b| ) we obtain the following representation for Q(u',u): 



Q(u',u) = ^ ( E KE n )S(E n+l -E n -u')S(E n+m -E n -u;)c n!n+l \ = AQ(0,w». (5.7) 



We now apply the Kirkwood approximation [23| to the correlation function, Eq. (|5.5b ), 
in the following form: 



„, n , , {p(0,r)p(u',r)){ P (0,r)p(u)){ P (u;',r)p(uy , 

Q(0, J y u; r) = A ^ ^ '- + B(0,u',uj; r) 5.8 

P 3 

where the numerator is chosen to incorporate all pairwise correlations, and the denominator 
ensures the correct limiting value as each energy difference tend to infinity. The correction B 
is small if the approximation is good, which is expected to be the case unless E", E' = E"-Ho' 
and E = E" + oj are all close together. Now from the definitions of p(E) and p(E,r), Eqs. 
fl2.2Q and (|5.6| ), the averages in Eq. (|5.8|) can be expressed via TLCF, R(u), 



[p(E + u , r) p(E)) = {p(E + w) p{E)) = p 2 (1 + R(u)) , 
and the Fourier transform of the return probability, p(u), 

Jd d r (p(E + u, r) p(E, r)) = \p(u) + 2>xhp\ = ^c(uj) 



16 



where we have taken into account that p(u) is expressed via irreducible part of the LDoS 
correlation function, and used the relation fl3.10p . 

Thus the Kirkwood approximation of Eq. ( |5.8| ) yields the expression 

P 



Q(u',u) = ^{c(u')[l + R(u)+R(u'-u)}} 



+ ^-c(u')R(u)R(uj'-uj) + A Jd d rB(0,u',u};r). 

The first line of the above equation represents the terms included in the Brownian motion 
formula, and the second line is the correction. We can now rewrite the exact relation Eq. 
Q4.7| ) using Eq. ( |5.4| ) and the factorization given by the first line of Eq. ( |5.9| ) to obtain 

m " fej -JdFcIF) (5 ' 10) 



where the remainder A(t) is given by the Fourier transform in the form ( |5.4a| ) of the second 



line of Eq. (|5.9| ). Note that in deriving Eq. ( |5.10|) , we, as above, have neglected a non 



integrated S(t) function which does not contribute to the final result. 

If we ignore the remainder A(t), Eq. ( |5.10|) is equivalent to Eq. (|5.3|) and thus to the 



Brownian motion result, Eq. (|2.5| ) We see that the use of the Kirkwood approximation is 
just a more formal way of doing the factorization discussed previously. The factor of 2 
that occurs in Eq. ( |5.2|) emerges naturally, as the two equivalent two-level correlations are 
automatically taken into account. It is obvious that the second line of Eq. (|5.9| ) could be 
small only if \u\, \u'\, \u — uj'\ A. Thus we could expect A(t) to be small, and Eq. (|2.5|) 
to be valid only outside of the quantum regime, i.e. for t < t H . 

What we would like to be able to do is to find out how big the remainder term A(t) 
actually is in this non-quantum limit. We can see that A(t) consists of two parts:- the first 
involves a product of three two-level correlation functions, and is the Kirkwood approxi- 
mation's attempt to represent three-level correlations in terms of two-level correlations; the 
second is the correction to the approximation itself. To proceed we will introduce diagram- 
matic perturbation theory in the next section, rewrite our factorization in this language, 
and hence discover what is left over. 



VI. DIAGRAMMATIC ANALYSIS 

Our aim now is to check the validity of Eq. ( |2.5| ) beyond the trivial diagonal approxi- 
mation in which one can just neglect the i-dependence of the denominator, and substitute 
the classical limit of p(t) into the numerator of this expression. The simplest way of doing 
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this is to calculate both K(t) and p(t) to higher order in perturbation theory, and to com- 
pare directly the r.h.s. and l.h.s. of Eq. ( |2.5| ). It is more instructive, however, to analyze 
diagrammatically the exact relation (|4.7| ) and the decoupling procedure based on Eqs. ( p.4| ) 
and ( |5.9| ). In this way we will not only check the accuracy of Eq. ( |2.5| ), but arrive at an 
alternative factorization given by Eq. ( |2.9| ). 

Let us note that the role of higher order corrections is different in d = 2 and d > 2. In 
the 2d case they are universal in the sense that they are due to diffusive motion of electrons 
throughout the whole sample and almost insensitive to details of motion at the ballistic scale. 
In d > 2, corrections are mainly due to the motion at the ballistic scale, and proportional 
to powers of an additional small parameter {t/t e \) d l 2 ~ l (as well as to powers of the standard 
weak disorder parameter - inverse dimensionless conductance). These corrections are not 
only small but of no particular interest as they do not drive the system from weak to 
strong-disorder regime. In contrast to this, the 2d corrections do describe crossover from 
the weak to strong disorder, and are widely believed to be more relevant for the vicinity of 
the metal-insulator transition for d > 2 than those calculated in the metallic limit directly 
in d > 2 dimensions. Moreover, as in the diagonal approximation K(t) oc tp(t) ~ const for 
t 4rg at d = 2, the TLCF vanishes in this approximation in the diffusive regime (E E c ). 
Therefore, in this regime the first non- vanishing higher-order contribution governs the main 
effect rather than describing some correction. For all these reasons, we will consider mainly 
the 2d case in this section. 



A. General Relations 

To be able to work in complete generality we will rewrite the exact equation ( f4.7|) in terms 
of Green's functions that can then be expanded using standard diagrammatic methods. The 
form of the diagrams for K(t) is by now well known, but those for p(t) and the three-level 
correlator Q(t) are less familiar. 

We start with the standard expression for DoS in terms of exact Green's functions where, 
as well as elsewhere in this section, we use the units h — 1: 

p(E, r) = \G R (E; r, r) - G A (E; r, r)j , (6.1) 

where the retarded, G R and advanced G A Green's functions are defined by 

G R,A (E; r, r') = V ^^"^ . (6.2) 

For the density- density correlation function we get the formula 
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= (^) Jd d r Jd d r f (Q(0;r,r)Q(uj;r',r')), (6.3) 

where Q = G R — G A . The diagrams for this then consist of two loops, each with an external 
vertex (with coordinates r and r'), as shown in Fig. (4a). Averaging over the disorder 
ensemble then leads to the presence of impurity lines both within a loops and across loops. 
Only the connected diagrams contribute to R(uj), while the trivial unconnected diagrams 
are cancelled by —1 in the definition of R(u), Eq. ( |2.3|) . 
Similarly, p{uj) is given by 

PM = jd d r (p(E + u, r)p(E, r)> = ~ Jd d r (g(u; r, r)0(O; r, r)) . (6.4) 



The difference between the formulae ( B.3|) and (6^4) is that in the latter the densities of states 
are evaluated at the same point r in space rather than at different points. In diagrammatic 
terms this means that p(u) consists of two loops connected to a single external point r, as 
shown in Fig. (5). 

In the diffusive regime of the disordered metal it is impurity ladders describing diffusive 
motion which give the important energy dependent contributions. They must have both a 
R and A line - so only RA and AR diagrams can contribute. Since the latter are complex 
conjugates of each other it follows that we can write the expressions for R(u) and p(u) in 
the following form: 

R(w) = ^^e Jd d r Jd d r f (G R (Lu;r,r)G A (0;r',r')) , (6.5) 
p(u) = - sfte Jd d r (G r (uj; r, r)G A (0; r,r)) . (6.6) 

The impurity ladders (Fig. 4b) which connect the two loops describe either diffuson or 
Cooperon propagators given in momentum space by 

D( g;ti ;) = — L- 1 , . (6.7) 
2npt el Dq 1 — iuo 

We will consider both the cases with and without time-reversal symmetry, referring to them 
as the orthogonal case (j3 = 1) and unitary case (j3 = 2). Diagrammatically, the latter 
differs from the former by the absence of the Cooperon contributions which are time-reversed 
impurity ladders. 

The loops could be connected by arbitrary number of the ladders. In the lowest order 
there are the two contributions to R(u>) shown in Fig. 4a, and one contribution to p(u) (Fig. 
5). The dominant contribution to R(u) - which is called Altshuler-Shklovskii diagram - has 
two impurity ladders; the diagram with only one ladder gives much smaller contribution. 
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(d) 

FIG. 4. (a) The two one-loop diagrams for R{ui) in a disordered metal. The diagram on the left 
gives the (dominant) Altshuler-Shklovskii contribution; that on the right is smaller by a factor At, 
and is usually ignored, (b) Definition of the impurity ladders occuring in (a), (c) Rewriting the 
impurity ladder as an effective propagator, (d) Rewriting the diagrams in (a) using the notation 
of (c). The electron Green's function lines end up in the so-called Hikami boxes. 
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FIG. 5. The one-loop diagram for the quantum return probability p(u>) in the ladder and 
Hikami-box representations. 

Thus the perturbation order of a diagram is not determined by the number of ladders. 
A standard way [24,p5| to determine this order, and to make the calculation of diagrams 
more straightforward is rewriting them in the form where impurity ladders are represented 
as propagators (a wavy line in Fig. 4c); this is convenient since the ladders involve small 
momentum, gl < 1, and energy, ut e i <C 1. All other Green's function lines are absorbed 
into effective interactions between propagators (as shown in Figs. 4 and 5) known as Hikami 
boxes. Then the order of a diagram is the number of independent momenta occuring in 
the propagators, which is just the number of loops made of the wavy lines. (Each box 
corresponds to a single spatial point r and can be thought of as being contracted into a 
point; thus any diagram would consist of some number of wavy loops). The loop-expansion 
parameter here is 1/g where g is the dimensionless conductance of the sample. 

The one-loop diagrams for R(u) are shown in Fig. 4d. The two external vertices may 
be in separate (odd) boxes, or the same (even) box; we ignore the latter since they are 
smaller than the former by factor At e i ~ {^/L) d in the same order of perturbation theory. 
The remaining Hikami boxes that do not contain an external vertex have an even number 
of sides. 

Rewriting diagrams for p(u>) in terms of Hikami boxes we find that the external vertex 
(corresponding to the point r in Eq. ( |6.6[) ) has two boxes connected to it - these can either 
both be odd boxes, or both be "2-gons" (larger even boxes are not allowed since one can 
always string impurity ladders across them to recover the "2-gon" structure). All other boxes 
have no external vertex and so have an even number of sides. Note that all the Hikami boxes 
must be "dressed", i.e. single-impurity lines should connect (without mutual intersection) 
those non- adjacent sides of the boxes which have the same analyticity (both R, or both A) 
25] . In all diagrams below the boxes are assumed to be dressed. The one- loop contribution 
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(a) (b) ( C ) 




(d) (e) (f) 

FIG. 6. Two loop order diagrams to: R{uj) (a, b, c); p(tu) (d, e); and Q(u>) (f). All diagrams 
contribute in the orthogonal case, while only (d) and (f) contribute in the unitary case. 

to p(uj) is given in Fig. 5 in both the ladder and Hikami representations. All the two-loop 
contributions to p(w) and R(u) are given in the Hikami representation in Fig. 6. 

Before comparing the higher-order contributions made to Eq. fl2.5|) by the diagrams for 
R{ui) andp(u;), we consider the three-level correlator Q(uj',uj). Starting from Eq. (|5.7| ), and 
rewriting the densities of states via Eq. ( p.ip we get 

Q{<J,u) = a(J^) Jd d r jd d r'(Q{tyr,r)Q{u'-r,r)Qtuj-r'y)) (6.8) 

Following the procedure used for R(u) and p(u>) above we see that in the ladder representa- 
tion the diagrams for Q(u',u) consist of three loops, two of which are joined at the external 
vertex r, and the third having the external vertex r'. There are now several classes of dia- 
grams that are not fully connected, as shown in Fig. (7a). (The shaded strips there include 
symbolically all possible combinations of impurity ladders). The most trivial has all loops 
unconnected and yields a constant term. The two others are reducible as they have only two 
out of the three loops connected by impurity ladders, and contribute terms proportional to 
p 3 R(u>), p 3 R(u>' — uj) and pp(u'). These reducible contributions add up to give 

Quc(u', w) = V^A {2tt P [1 + R(lo) + R(u'-u)] + p(co')} (6.9) 
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We see that the r.h.s. of the above is similar , but not identical to, the first line of the r.h.s. 



of Eq. ( |5.9| ). The difference is that the is no p(o/) multiplying the R(cu) and R{uo' — to) in 
Eq. ( p.9[ ), only its constant part 2np. If we substitute this term into Eq. ( |4.7| ) we obtain 
K(t) = (2nh)^ 1 tp(t) . In other words, the reducible diagrams of the three- level correlator 
Q are enough to reproduce the diagonal approximation of semi-classics. In order to recover 
our formula, Eq. ( |2.5| ), we will need to analyze the irreducible contributions to Q given 
schematically in Fig. 7b. 

If we look at the irreducible contributions the first thing to note is that each diagram 
can only have ladders between two of the three pairs of loops. This is because ladders are 
always between an R and an A line, and so in an irreducible diagram either two loops are 
R and the third A or vice versa. There are therefore 6 possibilities which come in complex 
conjugate pairs: AAR and RRA; ARA and RAR; RAA and ARR. The irreducible part of 
Q(u)',u) can then be written in the form 

Qd^u) = ^-^p Jd d r |A'3m{(G A (0;r,r)G^( W ;r,r)G»>')) 
+ (G A (0; r, r)G R {u; r, r)G A {u'; r', r')) + (G*(0; r, r)G A (u; r, r)G A (u'; r' , r'))} 

When rewritten in the Hikami-box representation the diagrams for Q consist of the two boxes 
connected to external vertex r found in p(u) diagrams, plus a single odd box connected to 
external point r' as found in R(u) diagrams; these are then held together by wavy lines and 
even boxes with no external vertices. Qaar diagrams cannot have the "2-gon" structure in 
the part connected to external point r since such a structure ends in ladders, which cannot 
happen here as both lines of the ladder would be A. Qara and Qraa diagrams can have 
this "2-gon" structure. 

If we are to re-derive our main result Eq. ( |2.5|) diagrammatically we will need to show that 
the connected part of Q(u)', u) factors into a product of R and p. This can be seen by recalling 
the derivation of Eq. (|2.5|) via the Kirkwood approximation, where such factorization occurs 
in the first line of Eq. ( |5.9| ). The fact that only two pairs of loops can have ladders between 
them suggests a possible reason for this to occur. At this point we note that we are interested 
not in Q(uj', uS) itself but rather in Q(t) which is derived from it by Fourier transforms as in 
Eq. ( |5.4a| ). The analyticity properties of Q(u', u>) diagrams will lead to some types giving no 
contribution to Q(t). Eventually, a factorization emerges from this analysis which is different 
to that suggested in Eq. ( |5.9|) , but which, remarkably, yields the same spectral form factor 
in two dimensions, to third-loop order. 
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(b) 

FIG. 7. (a) The reducible contributions to the three-level correlator Q(to', to). The first diagram 
yields no contribution to the function Q(t), whilst the other two yield equal contributions to 
Q(t) that are proportional to K(t). (b) The irreducible contributions to the three-level correlator 
Q(uj',uj). For analyticity reasons the last two diagrams give no contribution to the function Q(t), 
so we need consider only diagrams of the first type. 
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B. Analytical properties of the three-level correlation function 



First let us prove that only the Qaar diagrams survive, starting with the formula for 
Q(t) given by Eq. 



(6.11) 



Q(t) = j'dt'f doo' f doje~ iJt ' e^ l Q{oj\oj) 



o 



doo' I duo — Q(oj', to) 



Taking the Fourier transform of Q(t) then yields 

D icut p i{iu—Lo')t 



J dte** J dJ I doo — Q(u',u) 

= 2tt I dJ fduj— \8(uJ + u) -5{uj + u-uj')]Q{uj',uj) (6.12) 
d 

—2m / — - [Q(oj', —uf) — Q(lu', oj' — uf)] 



Now consider the analyticity properties of the ARR, ARA and AAR respectively. They can 
be written in the form 

Qraa(^', lo) = /(-u/, -u) 

Qara(u', oj) = f(oo\ oj' - oj) (6.13) 
Qaar(^', oj) = f(oj, oj - oj') . 

where in each case f(ui,U2) is a function analytic in the u.h.p. for both arguments. For 
Qraa Eq. ( |6.12|) gives 



-2m f^- [f(-oj',uJ] - f(-oj', -oj'+Lj)] = -27T 2 [/(0, -oj) - /(0, -co)} = , (6.14) 

J OJ 

since upon closing each term in l.h.p. we only get contributions from the pole at oj' = 0, 
and these cancel. Similarly for Qara we get upon closing in the u.h.p. 

du' 



-2m I ^ [/(a/, oj' + cj)- f(oj', 57)] = 2vr 2 [/(0, 57) - /(0, 57)] = . (6.15) 

OJ 



Finally for Qaar we get 

-2m f^- [f(-ZJ, -oj' - oj) - f(u' - oj, -57)] = -(2n) 2 f(-ZJ, -oj) . (6.16) 

J OJ 

where we closed first term in l.h.p., second in u.h.p. to yield contributions that add up. 
The above term can be rewritten in the form — (2w) 2 Q aar{u' — 0, —57), so it follows that 
diagrams for Qaar can yield Fourier transform of Q(t) directly. We have therefore verified 
our previous assertion, and moreover have derived the contribution to the Fourier transform 



25 



of Q(t) coming from the Qaar diagrams. Diagrammatically the above means that all ladders 
in Qaar diagrams will have the same energy dependence. 

Let us discuss what the above means for our factorization hypothesis. The reason we 
expected that it would be the Qrar and Qarr diagrams that survived is that we can 
envisage a natural factorization into terms with energy dependence of the form R(u)p(u') 
and R{uj' — uj)p(uj'), which are exactly what is needed to reproduce Eq. (|2.5|). The fact that 
it is Qaar which survives means that our factorization hypothesis must be altered to have 
an energy dependence R(u)p(u), which leads to a convolution in t-space. More precisely, 
we expect the a;-space factorization to have the form 

Q {u) = J^h p{uj) R {uj) . (6 .17) 
Here we have introduced the analytical functions Q(oo), p(u>) and R(u) as follows: 

Q(u) = Qaar(0, w) ; 

p{u) = 2^ep(to); (6.18) 
R(u) ee 2^eR(u). 

Note that diagrammatically we calculate exactly these analytical functions. Before giving 
the results for higher-loop contributions, note that calculations of R(u) are considerably 
simplified due to the fact that it can be obtained as the second derivative of "free energy" 
F(w): 

Here F{u) is given by the sum of all diagrams that have no external vertices, and the 
coefficient of proportionality is chosen so as to make F(uS) dimensionless. 



C. Calculation of diagrams up to three-loop order 

As a result of the above discussion we can now lay out the diagrammatic program ahead. 
First, we will calculate the contributions to F(u), p(u), and Q{uS) up to three-loop order. 
Then we will first verify Eq. (|4.7|) (which has been obtained assuming spectral homogeneity), 
and then check the factorization scheme of Eq. fl6.17p which is dictated by the analytical 
structure described above. This factorization is equivalent to that given by Eq. ( |2.9| ). Finally, 
we will show how, having verified this factorization up to the three-loop order, to check the 
factorization that comes out of the Kirkwood approximation, Eq. ( p.9| ), and is equivalent 
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to the result of the Brownian-motion model, Eq. ( |2.5|) . We will show that up to the third- 
loop order both the factorizations are identical and exact. Therefore, Eq. ( |2.5| ) obtained 
within the Brownian motion picture turn out to be correct to quite a nontrivial order of 
perturbation theory. 

Let us introduce a notation that is convenient for further analysis: 

V x ... n = D( qi + . . . + q n f - iu (6.20) 

In particular, V\ = Dq\ — iu. 

Now, two- loop order contributions corresponding to the diagrams in Fig. 8(a) (equivalent 
to those in Fig. 6a-c) and in Fig. 6(d-f) may be written as 



2 



p 2 (u) = (l) *»e E j^V + 2(2 - 0)*±n±»L) ; (6.21) 



2\ 2 A ^ 1 



W = b Yd £ 



Here we note the following. The diagrammatic approach could be used both in the diffusive 
regime, uj ^> E c , where all the sums above should be replaced by integrals, and in the 
ergodic regime, E c ^> to 3> A, where only contribution of zero mode (all q = 0) survives. 
In the diffusive regime, regularization of divergent integrals is required. Although we do 
not explicitly calculate contributions of all diagrams below, in all algebraic manipulations 
we use dimensional regularization near d = 2. These manipulations involve dealing with 
large q ~ and so all diagrammatic expressions we give here are not directly valid in the 
ergodic regime. However, it is straightforward to verify the accuracy of our factorization 
scheme in the zero-mode regime as well. 

The three-loop contributions to F(u), p(uj), and Q{uS) are shown in Figs. 8, 9, and 
10, respectively. In Fig. 8a we have also drawn the two-loop contribution to F(u) which 
is equivalent to the three two-loop contributions to R(u) shown in Fig. 6. Note that in 
the three-loop order there are 41 diagrams that contribute directly to R{uj), so that to 
have instead only the 5 diagrams for F(u), as in Fig. 8, is a considerable simplification. 
Nevertheless, these three-loop results are quite bulky, and we list contributions of each 
diagram in the three tables. 
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(2) (3a) (3b) 




(3c) (3d) (3e) 

FIG. 8. Two and three-loop order diagrams for the function F(u). The two-loop diagram is 
equivalent to those in Fig. 6(a-c), and contributes only to the orthogonal case. All 5 three-loop 
diagrams contribute to the orthogonal case, while only 3a and 3b contribute to the unitary case. 




(a) (b) ( C ) (d) (e) 




(f) (g) (h) (i) (j) 



FIG. 9. Three-loop order diagrams for the quantum return probability p(u>). All 10 diagrams 
contribute in the orthogonal case, while only the first 2 contribute in the unitary case. 
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(a) (b) (c) (d) 




CD (k) (1) (m) 



FIG. 10. Three loop order diagrams for the function Q{ui). All 13 diagrams contribute in the 
orthogonal case, whilst only the first 4 contribute in the unitary case. 
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3F a -F d --2 ^ VMP* 

91,92,93 ' 1' 2' 3 

„ i ^ {V 12 + V 23 ) 2 



4 91^93 nV 2 V 3 V 12 3 



91,92,93 



1 (^12+^23)(^13+^3 2 : 



2 91,92,93 V{P 2 VzVyf 

Table 1. Three-loop order contributions to F(u>): an overall factor of (A/27r) 2 should 
be attached to each contribution. 



1 ^ 2(V 1 + V 2 + V 3 )+3iu) 



Pa = ~Pe=Pf = -2 E r- 

Pt> = E 



2 1 ^ ViVoV 3 

91,92,93 ' 1 ' 2' 3 

gig + ^ 23 ) 2 
, 3 V\V 2 V:{P\x; 



, v (Pi+P 2 + ^)(Pi + P 3 +^) 

91,92,93 ' 1 ' 2' 3 

, v (Pl+P 2 + ^)(^2+^3 + H 
^ d — 4 2^ ^2^2 



p 9 = z E 

9i,92,' 

Pfc = 4 E 



-p2-p2-p 

91,92,93 1 2 ' 3 

QPl2 + 7 7 23)(7 7 r3 + P 3 2) 
f/ 1 ,(/2 ■('/:■; ^1 '2^3' 123 

(Vx +V 2 + iuj) 



91,92,93 ' 1 ' 2' 3 

91,92,93 ' 1' 2' 3 
gl2 + ^23) 

V{P 2 V 3 V l23 



pj= E 



91,92,93 



Table 2. Three- loop order contributions to p(oj): an overall factor of L d (A/27r) 2 
should be attached to each contribution. 
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Qa — %Qb — Qh — Qi — Qj — Qk — 2i £ -p2-p -p 

91,92,33 ' 1 ' 2/ 3 



Qc = Ql=i 

(Pl2 + P 23 ) 



Qd — Q 

Qe 

Qf 

Qg 



V{P 2 V Z V 
-2i e 



91,92,93 ^1 ^2^3^123 

(Pi + P 2 + zw) 



4z V 

91,92,93 1 2 ' 3 

(Pi + P 3 + iw) 



-2i £ 



91,92,93 



(Pi + V 2 + iu) 



91S93 PiP2 2 P3 



Table 3. Three-loop order contributions to Q{u): an overall factor of (A 2 /27rL d ) 
should be attached to each contribution. 



The labels in the tables correspond to those in Figs. 8-10, and in the figure captions 
we describe which diagrams are made of diffusons only and thus contribute in the unitary 
{(3 = 1) case. Obviously, all the diagrams contribute in the orthogonal (/3 = 2) case. Results 
for the symplectic symmetry class {(3 = 4) are practically the same as for the orthogonal 
case but we do not list them here to avoid complications with the coefficients. 

The starting point for our diagrammatic analysis is Eq. (|4.7|) which we now rewrite as 
follows. The Fourier transform, Eq. (|5.4j) , of the function Q(t) is split into the reducible and 
irreducible parts (Eqs. ( |6.9|) and ( |6.10| ), respectively). As only Q(w), Eq. (|6.18|) , contributes 
to the irreducible part, we obtain 







d(iu>) 



L d d(iu 



F(u)-p(u) 



-Q{u) (6.22) 

7T 



This is the exact relation which holds to all orders in perturbation theory. It allows us to 
check the accuracy of our diagrammatics up to three loop order for both the unitary and 
orthogonal cases. We do this by first substituting the two-loop results of Eq. (|6.21 ) - which 



is quite straightforward, and then the three-loop data from the tables which requires some 
significant algebra. We verify that this identity holds with our diagrammatic results which 
gives us confidence in their accuracy. 

The next step is to see whether our projected c<j-space factorization occurs to two and 
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three loop order in both the orthogonal and unitary cases. To check Eq. ( |6.17|) up to two- 
loop order, we calculate both factors in the r.h.s. to the first order only. This is simple and 
yields 

so that the factorization ( |6.17D is exact in this order. The reason why this is so simple 
is that with only two momenta q\ and q 2 there is no way for the momenta to become 
"entangled", and so there is only really one possible functional form. The fact that the 
numerical coefficients match up exactly is the important thing. The three loop case is more 
involved because now the entanglement can occur, and this leads to the factorization not 
being exact. For both the unitary and orthogonal case, however, we get the same functional 
form in the remainder, 

i(2n) 2 r ~ ~ i 

Q 3 (u) + + R 2 (uj)Pi(uj) 

2zA 2 f 4 2(V 12 + V 23 ) 1 1 



*l3 2 L d qi ^ q3 { (TJ 2 V 2 V 3 (n) 2 V 2 V 3 V 123 Vi ^3 P123 J ' 
The remainder can then be algebraically manipulated to give 

8 (6.23, 



37T/W d(iu) q ^ m \V x V 2 V 3 V 123 

At the three loop level we see that the factorization is not exact, but that the remainder term 
is simple and its momentum integrals are fully convergent. Certainly for the orthogonal case 
many terms have been removed to yield this remainder. We next note that in the special 
case of two dimensions this remainder is zero, and the factorization is exact. This is because 
for d = 2 dimensional analysis shows that the term inside the derivative is a constant - it is 
of the form (— and no logarithmic singularities are present - and so one gets zero upon 
taking derivative. Even if the remainder were not able to be written as a derivative, but 
just as a sum of terms with no logarithmic singularities, it would yield zero. This is because 
dimensional analysis shows that the result is of the form a/ (—iu), where a is a real constant. 
Since we have to take the real part of this to get R(u) this would give zero. However in this 
case there would be a constant contribution to the Fourier transform K(t), because a/ (—iu) 
does actually have a real part proportional to 5 (id) - this is exactly what happens in the 
case of the 1-loop contribution in d = 2. We note that our factorization is exact to 3-loop 
order in d = 2 even in the sense of getting the constant term in K(t) correct. 

It seems to us that the exactness of factorization up to three loop order in d — 2 for 
both orthogonal and unitary cases is no accident, and we conjecture that this result persists 
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to all orders in perturbation theory. Obviously such a conjecture cannot be proved using 
order-by-order analysis (although, of course, it could be disproved this way), so any attempt 
to verify this will require analysis of the structure of R(u), p(u), and Q{uS) diagrams. 

D. Comparison of the Two Factorization Schemes 

In this section we will compare the two factorization schemes that we have introduced in 
this paper: the t-space scheme that arose from consideration of the Brownian motion model 
of section (III), and the u;-space scheme that arose in the diagrammatic analysis of section 
(VI). We have shown that the cj-space scheme is exact to 2-loop order in all dimensions, 
and to 3-loop order in the 2d case. We will now examine the validity of the t-space scheme. 
This involves very little extra work because most of the algebraic manipulation has already 
been performed in the cu-space analysis. We first compare the two relations by writing both 
of them in t-space to yield 

K(t) - (2irhp)-Hp(t) = -(nhpy 1 f dt'K(t')p(t - t') (6.24a) 

K{t) - [2ithp)-Hp{t) = -(nhp^Kit) f pit') (6.24b) 

Jo 

where, of course, we know the region of validity of the first formula. To investigate the 
t-space factorization we need only compare the r.h.s. of the above equations. In the 2- 
loop case the K(t) and p(t) in the r.h.s. will both be of 1-loop order, and we know that 
Pi{t) oc t~ d l 2 and Ki{t) oc t l ~ d l 2 . We find that the two equations above only agree for d = 2, 
where K~i{t) = 1. The t-space factorization is therefore exact to 2-loops only in 2d, and we 
restrict ourselves to the 2d case from now on. 

To look at the t-space scheme to 3-loop order in 2d we note that we can have the 
combinations K ± , p 2 , and K 2 , p\. For the unitary case K 2 (t) = 0, so this leaves only the 
first contribution, and since Ki(t) = 1 the two factorizations become the same. For the 
orthogonal case we have to look at the second contribution. We find that the two equations 
differ by a constant. For the purpose of calculating the asymptotic behaviour of R(cu), 
constant terms in K(t) have no effect, so that the t-space factorization works up to 3-loop 
order in 2d for both orthogonal and unitary cases. 

At this point it seems that the tu-space scheme may be perturbatively slightly more 
accurate than the t-space scheme in that it is correct in 2-loops for all dimensions, and at 
3-loop order it gets the constant term in K(t) correct in the orthogonal case. However we 
are still justified in saying that both schemes are correct to 3-loop order in 2d. 
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VII. SUMMARY 



In this paper we have examined spectral correlations in disordered conductors, starting 
from the idea (Eq. |4.2|) that two samples with impurity configurations differing by an in- 
finitesimal amount should be statistically equivalent. In the first instance, this equivalence 
yields an identity relating two-point to three-point correlation functions; it is useful only if 
one can decouple the latter. We argue that a decoupling based on the Kirkwood superposi- 
tion approximation is physically reasonable provided one is interested only in correlations at 
scales large compared to the mean level spacing. Within the approximation, we express (Eq. 
|2.5| ) both the non-parametric and the parametric spectral form factor in terms of the quan- 
tum return probability for a spreading wavepacket. We test the decoupling by calculating 
corrections, using the standard diagrammatic perturbation theory for disordered conductors 
to expand about the metallic limit in inverse powers of the dimensionless conductance, g. 
We show that in two-dimensional systems, the case of greatest interest, there are no cor- 
rections to order g~ 3 . We believe that the results we obtain from this approach should be 
useful rather generally, and especially when a diagrammatic analysis is not straightforward, 
as at the Anderson transition; the implications of our work in that regime will be discussed 
elsewhere. 



ACKNOWLEDGMENTS 

We thank V. E. Kravtsov and B. D. Simons for numerous useful discussions. Support by 
the EPSRC under grants Nos. GR/GO 2727 (J.T.C.) and GR/J35238 (I.V.L. and R.A.S.) 
is gratefully acknowledged. I.V.L. acknowledges the hospitality of the ITP in Santa Barbara 
where part of this work was performed, and partial support by the NSF under grant No. 
PHY94-07194. 



34 



REFERENCES 

[1] E. P. Wigner, Proc. Cambridge Philos. Soc. 47, 790 (1951). 
[2] F. J. Dyson, J. Math. Phys. 3, 140 (1962). 

[3] M. L. Mehta, Random matrices (Academic Press, Boston, 1991). 

[4] L. P. Gor'kov and G. M. Eliashberg, Zh. Eksp. Teor. Fiz. 48, 1407 (1965) [Sov. Phys. 
JETP 21, 940 (1965)]. 

[5] K. B. Efetov, Adv.Phys. 32, 53 (1983). 

[6] B. L. Altshuler and B. I. Shklovskii, Zh. Eksp. Teor. Fiz. 91, 220 (1986) [Sov. Phys. 
JETP 64, 127 (1986)]. 

[7] N. Argaman, Y. Imry, and U. Smilansky, Phys. Rev. B 47, 4440 (1993). 

[8] J. T. Chalker, I. V. Lerner, and R. A. Smith, submitted to Phys. Rev. Lett., 1996. 

[9] F. J. Dyson, J. Math. Phys. 3, 1191 (1962). 

[10] P. Pechukas, Phys. Rev. Lett. 51, 943 (1983). 

[11] M. Wilkinson, J. Phys. A 21, 4021 (1990), P. Gaspard et al, Phys. Rev. A 42, 4015 
(1990), J. Goldberg et al, Nonlinearity 4, 1 (1991). 

[12] A. Szafer and B. L. Altshuler, Phys. Rev. Lett. 70, 587 (1993). 

[13] B. D. Simons and B. L. Altshuler, Phys. Rev. Lett. 70, 4063 (1993). 

[14] B. D. Simons and B. L. Altshuler, Phys. Rev. B 48, 5422 (1993). 

[15] T. Yukawa, Phys. Rev. Lett. 54, 1883 (1985). 

[16] C. W. J. Beenakker, Phys. Rev. Lett. 70, 4126 (1993). 

[17] A. I. Larkin and D. E. Khmelnitskii, Uspekhi Fizicheskikh Nauk 136, 536 (1982). 

[18] D. E. Khmelnitskii, Physica B&C 126, 235 (1984). 

[19] V. E. Kravtsov and I. V. Lerner, Phys. Rev. Lett. 74, 2563 (1995). 

[20] In the ergodic regime, this expression gives only the envelope of the TLCF as the 
oscillations in R(s) || are due to singularities in K(t) near t H , which cannot be accounted 
for semi-classically. 

[21] F. J. Dyson, J. Math. Phys. 13, 90 (1972). 



35 



[22] L. A. Pastur, Th. Math. Phys. 10, 67 (1972). 
[23] J. G. Kirkwood, J. Chem. Phys. 3, 300 (1935). 

[24] L. P. Gor'kov, A. I. Larkin, and D. E. Khmel'nitskii, Pis'ma v ZhETF 30, 248 (1979) 
[JETP Letters 30, 228 (1979)]. 

[25] S. Hikami, Phys. Rev. B 24, 2671 (1981). 



36 



